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Abstract 

One way to predict hurricane numbers would be to predict sea surface temperature, and then pre- 
dict hurricane numbers as a function of the predicted sea surface temperature. For certain parametric 
models for sea surface temperature and the relationship between sea surface temperature and hurri- 
cane numbers, closed- form solutions exist for the mean and the variance of the number of predicted 
hurricanes, and for the standard error on the mean. We derive a number of such expressions. 

1 Introduction 

One way to try and predict future hurricane numbers is to predict sea surface temperatures (SST), and 
then to predict hurricane numbers as a function of the predicted SSTs. If both the prediction of SSTs 
and the model that relates SSTs to hurricane numbers are probabilistic then the resulting probabilistic 
prediction of hurricane numbers can, in general, only be derived using numerical methods. However, in 
certain cases a lot of information about the predicted distribution of hurricane numbers can be derived 
analytically. This includes estimates for the mean number of hurricanes, the variance of the number of 
hurricanes, the standard error on the estimate of the mean number of hurricanes and the linear sensitivity 
of the mean number of hurricanes to changes in the mean and variance of SST. 
In this article we derive a number of such relations, for the following two cases: 

1. we predict the SST distribution, and then use these predicted SSTs to predict either basin hurricane 
numbers or landfalling hurricane numbers. 

2. we predict the SST distribution, use this to predict the distribution of basin hurricane numbers, 
and then use a further relationship to predict landfalling hurricane numbers from basin hurricane 
numbers. 

The assumptions we make to render this problem tractable are as follows: 

• the sea surface temperature is taken as normally distributed (although in several of our derivations 
we relax this assumption and consider a completely general distribution of SST with known moan 
and variance) 

• the relationship between sea surface temperature and hurricane numbers is taken as either (a) 
linear and normally distributed, (b) linear and poisson distributed, or (c) exponential and poisson 
distributed. 

• the relationship between hurricane numbers in the basin and hurricane numbers at landfall is taken 
to be linear and poisson distributed 

This article proceeds as follows. In section[5]we discuss our choice of models and assumptions. In section[3] 
we discuss the types of statistical models we will use, and the terminology we will use to describe them. 
In section 2] we describe our notation and what we need from the SST forecasts. 

In section [5] we derive expressions for aspects of the predicted distribution of hurricane numbers in the 
case where the relationship between SST and hurricane numbers is linear and normally distributed. In 
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section [5] we derive expressions for aspects of the predicted distribution of hurricane numbers in the case 
where the relationship between SST and hurricane numbers is linear and poisson distributed. In section[7] 
we derive expressions for aspects of the predicted distribution of hurricane numbers in the case where the 
relationship between SST and hurricane numbers is exponential and poisson distributed. 
In section[8]we derive expressions for aspects of the predicted distribution of landfalling hurricane numbers 
when predicted as a function of the number of basin hurricane numbers, in the case where the relationship 
between the two is linear and poisson distributed. 

Finally in section [9] we discuss how to predict landfalling hurricane numbers from basin numbers which 
are themselves predicted from SST by combining the relationships derived in sections [5l [5] and [7] with the 
relationships derived in section [HJ 



2 Comments on our choice of models 



This paper derives various mathematical expressions related to the prediction of hurricane numbers. It is 
not, however, a discussion of which models are actually appropriate to use to predict hurricane numbers. 
This latter question is a questio n we discuss at great l ength else where: for a di s cussio n of what models 
can be used to predict SST see Meagher and Jewsonl ( 2006f ) and Laepple et al. ( 20061); for a d i scussio n 
of what models can be used to relate SST to landfalling hurricane numbers see Binter et al.l ( 2006b % ) 



for a discussion of what models can be used to relate SST to basin hurricane numbers see Binter et al 



(|2006q ): and for a disc ussion of wha t model s can be used to relate basin hurricane numbers to landfalling 



hurricane numbers see lBinter et alJ (|2006al ). 

However, we now give a brief summary of the rationale behind our model choices, based on the results 
from these studies: 



we consider a normal distribution for SST because that seems to be the simplest reasonable model 

we consider a linear-normal relationship between SST and hurricane numbers because this is the 
simplest case analytically, even though the use of the normal distribution for hurricane numbers 
may not be reasonable if the number of hurricanes is small 

we consider a linear-poisson relationship between SST and hurricane numbers because this is the 
simplest case that has a reasonable distribution for hurricane numbers even in the situation in which 
the number of hurricanes is small. We are not concerned that use of a linear relationship could in 
principle lead to a negative value for the poisson parameter since this does not happen in practice 
with the data we are using. 

we consider an exponential-poisso n relationship between SST and hurricane numbers because this 
has been used previously (e.g. bv lElsner and Schmertman n (1993)), and is the stand ard way that 



statisti cians tend to use poisson regression. However, from the analysis we describe in lBinter et al 



(|2006bf ) we conclude that it is not possible to tell from the data whether this model is better or 
worse than the linear-poisson model. 

we consider a linear-poisson relationship between basin hurricane numbers an d landfalling hurrican e 
numbers because it is simple, it does well in our own tests with real data (see Binter et al. (|2006al )) 
and it includes the simpler model that consists of just a constant proportion. 



3 Comments on notation and statistical models 

Consider trying to build a statistical model that models some variables j/j as a function of some other 
variables Xi- One obvious place to start is standard linear regression, which can be written as: 

Ui = a + /3xi + ei (1) 

where the are taken to be IID and normally distributed with mean zero. From now on we drop the 

subscripts i for simplicity. 

We can also write this model as: 

y~N(a + px,o- 2 ) (2) 
In this paper we will call this the linear-normal model. 

This model can be generalised in a number of ways. One way would be to change the distribution from 
normal to something else. We could do this either by specifying a different distribution for the noise 



forcing e, or by specifying a different distribution for the response y. In general these are not equivalent: 
for instance if we were to specify that the noise should be poisson distributed then y would have a poisson 
distribution shifted by a + fix, and, conversely, if we were to specify that y should be poisson distributed 
then the noise would have a poisson distribution shifted by —a — fix. In our case we think we have more 
idea about the distribution of y than we do about the distribution of the noise: since y is modelling 
hurricane numbers we think that it will be close to poisson. So the first generalisation we consider will 
be to replace equation [2] with: 

y ~ Po(rate = a + fix) (3) 

In other words: we model y as poisson distributed with mean a + fix. We will call this the linear-poisson 
model. 

We can also write this model as 

y = a + fix + e (4) 

just as before, but we note that although it is still the case that E(e) — 0, the distribution of e is now 
somewhat odd, since it is a poisson distribution but shifted to have zero mean. As a result it is not 
common to write this model in this way (although there is nothing wrong with it). 
We could also write this model using conditional expectations: 

E(y\x) = a + fix (5) 

with the additional information that y is poisson distributed. 

We note that the linear poisson model does not often appear in statistics textbooks because for certain 
values of a, fi and x the rate of the poisson could be negative, which is impossible. We, however, take a 
practical approach: for our data this is not a problem, and what happens far outside the domain covered 
by our data is not relevant to us and can be ignored. 

Having discussed how to change the distribution used in standard linear regression, the second obvious 
generalisation would be to replace the linear function with something non-linear, and a common way to 
do this is to replace a + fix with exp(a + fix). For those concerned about the problem with negative 
poisson rates described above, this solves that problem. We can write this new model as: 

y ~ Po(rate = exp(a 4- fix)) (6) 

or 

y = exp(a + fix) + e (7) 

or 

E(y\x) = exp(a + fix) (8) 

where y is poisson distributed. 

In fact, statisticians often write this model in yet another way, as: 

\ogE{y\x) = a + fix (9) 

and call it the log-linear poisson regression model, although we find this nomenclature rather unhelpful, 
and we prefer to call this model an exponential-poisson regression model. 

One final comment is that one can discuss all the models given above in the context of a general class of 
models known as generalised linear models (GLMs), which cover almost any possible distribution for y 
and any possible non-linear monotonic function of a + fix. 

4 Setup and basic relations 

We will use the following notation: 

• s is the SST 

• s is taken as normal, with mean /i s and sd a s 

• n is the number of hurricanes 



• p() is used for probability densities 



• E() is used for expectations 

• when we need to distinguish between the number of basin hurricanes and the number of landfalling 
hurricanes we write these as rib and n;. 

• the relationships we use to model hurricane numbers as a function of SST give us the parameters 

a,(3,a 

• the relationships we use to model landfalling hurricane numbers as a function of basin hurricane 
numbers give us the parameters a',p', a' 

• the mean and variance of the number of hurricanes is written as /j,h and o 2 , with and o 2 for 
basin hurricanes and m and of for landfalling 

4.1 Calculating Means 

We will calculate the mean number of hurricanes using: 

Hh = E(n) (10) 

oo 

= ^2 n p( n ) (ii) 

°° roc 

= J2 n P{n\s)p{s)ds (12) 

n=0 J -°° 

/oo / 00 \ 
[Y,np{n\s)\ P {s)ds (13) 
-°° \n=0 / 



/oo 
E(n\s)p(s)ds (14) 
-oo 



4.2 Variances 

We will calculate the variance of the number of hurricanes using: 



o\ = E[(n-^ h f] (15) 

OO 

= ]T(n-^) 2 p(n) (16) 

n=0 

= Y,( n -^h) 2 / p(n\s)p(s)ds (17) 

/oo / oc \ 
PT^-^pHs) Ms)ds (18) 
-oo \n=0 / 

/oo 
^(n-^) 2 ^^ (19) 
-oo 

In the non-linear cases we will calculate the variances using: 

sigmal = E[(n-^ h f] (20) 
= E{n>)-ti (21) 

where 

oo 

E(n 2 ) = 5> 2 p(n) (22) 

°° /'OO 

= 5Z n / PH s M s ) ds (23) 



n=0 



n=0 



/ \^2n 2 p(n\s)\ p(s)ds 
Jn=o \ n=0 ) 



3C 



(24) 

E{n 2 \s)p(s)ds (25) 



n=0 

4.3 SST forecasts 

We assume that we have some method for the prediction of SST that gives us a normal distribution 
N(/j s ,a 2 ). The parameters and a s will typically have some uncertainty associated with them, which 
will depend on the prediction method being used. Some of the expressions we derive below need estimates 
of this uncertainty in order to derive estimates of the standard errors on our hurricane number predictions. 

5 The linear-normal model for the relationship between SST 
and hurricane numbers 

We now describe the first of our models for the relationship between SST and hurricane numbers, which 
is the linear-normal model. In this model we postulate a linear relation between SST s and hurricane 
numbers n, where n could be either the number of hurricanes in the basin, or the number at landfall. We 
assume that the distribution of hurricane numbers is normal. As discussed above, we include this model 
because it is perhaps the simplest model one might consider. However, the assumption of normality 
probably doesn't hold very well when the number of hurricanes is small: this is resolved in section [5] by 
replacing the normal distribution with the poisson. 
The linear-normal model can be written as: 

n = a + Ps + e (26) 
e ~ N^cr 2 ) (27) 

5.1 The predicted mean 

Taking expectations of equation [26] (over all realisations of e and s) gives us a simple expression for the 
mean number of predicted hurricanes in this model: 

Hh = E(n) =a + pE(s) =a + (3fi s (28) 

We see that the mean number of hurricanes is a linear function of the mean SST, and doesn't depend on 
the variance of SST a 2 . 

5.2 The predicted variance 

We can also derive an expression for the variance of the number of hurricanes a\ fairly easily. 
Combining equations |2T)1 and |2"51 we see that: 

n-n h = f3(s-n s )+e (29) 
(n-n h ) 2 = fJ 2 ( S -n s ) 2 + 2p(s-n s )e + e 2 (30) 
E((n~fi h ) 2 \ S ) = p 2 ( s -^ s ) 2 +a 2 (31) 

where we take expecatations over all realisations of e, for fixed s. 
This gives: 



a 2 = I E((n- n h ) 2 \s)p(s)ds (32) 
(P 2 (s- ( i s ) 2 +a 2 )p(s)ds (33) 
P 2 f(s- Li s ) 2 p(s)ds + ( <r 2 p{s)ds (34) 



H 2 a 2 + a 2 (35) 



This expression is easy to understand: the variance in the number of predicted hurricanes comes both 
from the variance in the SST prediction (scaled by /?) and the variance around the relationship between 
mean SST and the mean number of hurricanes. The variance of the predicted number (b) the standard 
error on the SST prediction (which should be given by the SST prediction routine). 



5.3 Standard errors 

We now derive an approximate expression for the standard error on the prediction of the expected number 
of hurricanes. 
We already have: 

Hh = a + P/j, s (36) 

Now consider small errors in the parameters a and (5 and a small error in the prediction of the mean 
SST. We can understand what errors this leads to in our prediction of expected hurricane numbers just 
by linearising: 



dm = da + + f3d^ (37) 
Squaring this and taking expectations gives: 

var(/x/j) = var(a s ) + ^i 2 var(/3) + /3 2 var(/x s ) + 2^ s cov(a, /3) (38) 

This gives us an approximation to the standard error on the expected number of hurricanes in terms of 
the standard errors on the parameters from the regression (which are given by most regression routines), 
and the standard error on the SST prediction (which should be given by the SST prediction routine). 



5.3.1 Higher order terms 

By using Taylor expansions, we can derive higher order terms in the expression for the standard error. 
Consider the function fih = f(a, /3, /i s ). Expanding this function in a Taylor series gives: 



d^h = -^—da + -—df3+- — dy, s (39) 
da op o/j, s 

+ 2W da + 2F* + 25^ (40) 

(41) 

= da + ii s dj3 + (3dfi s (42) 

+0da 2 + Od/3 2 + Qdul (43) 

+0dad/3 + ld/3dfi s + Od^ s da + (44) 

= da + fi s dp + (3d^i s + d(3d^i s (45) 

This is exact, since all subsequent terms are zero. Squaring and taking expectations now gives: 

var(d^/ l ) = var(a) + fi 2 YBx{(3) + 2^i s cov(a, (3) + /3 2 var(^i s ) (46) 
and we find that our original expression is actually exact. 



5.4 Linear sensitivity 

We can also consider the linear sensitivity of our forecast to changes in the mean and the standard 
deviation of our SST prediction. This can be useful to understand how a change in the forecast will 
create a change in the hurricane prediction. 
Again starting with 

[i h = a + Pn s (47) 

if we differentiate wrt [i s we get: 



(48) 



and 

^ = (49) 

0<7 S 

and we again see that the predicted mean number of hurricanes is independent of the variance in the SST 
forecast. 



5.5 Summary 

We now summarise the relations we have derived for the linear-normal model: 



fih = a + (3fi s (50) 
a\ = (3 2 a 2 s +a 2 (51) 
var(/ih) = var(a) + /J 2 var(/3) + /3 2 var(/j s ) + 2/j s cov(a,,S) (52) 

d t, - " (53 > 

We didn't actually use the fact that the SST was normally distributed to derive these relations: all we 
used were the values for the first two moments of the SST distribution. So all the above results hold for 
any SST distribution, given the first two moments. 



5.6 Alternative representation 

In practice, equation [52] is difficult to evaluate because the right hand side contains two large positive 
terms (var(a) and /i 2 var(/3)) and one large negative term (2/x s cov(a, /?)). Rounding error can easily cause 
the result to be negative (when it should be positive), or at least very inaccurate. 
We can avoid this problem by rewriting the original regression equation as: 

n = a + f3(s - s) + e (55) 

where s is the mean of the observed values of historical SST. Compared with the original formulation the 
value of (3 is the same but the value of a is now different. 
This then gives: 

Hh = a + (3{n s - s) (56) 

The expression for the variance, which doesn't depend on a, doesn't change. 
The derivation for the standard errors is based on: 

d/j,h = da + (fi s — s)d/3 + ftd(X s (57) 

giving: 

var(fih) — var(a) + (^ s — s) 2 var(/3) + /3 2 var(^ s ) (58) 
The cov(a, 0) term drops out because it is zero in this alternative representation. 



6 The linear-poisson model for the relationship between SST 
and hurricane numbers 

We now derive relations for the second of our models for the relationship between SST and hurricane 
numbers. This model is the same as the previous model, except that we replace the normal distribution 
with a poisson distribution. Because the expression for the mean number of hurricanes given the SST is 
still linear the results are rather similar. 
We write this model as: 

n = a + fis + e (59) 
ri ~ Po(rate = a + (3s) (60) 



Note that the variance of n given s is given by v(s) = a + (3s whereas in the previous model the variance 
was constant. 



6.1 The predicted mean 

Just as before: 



Hh = E{n) = a + f3p, s (61) 

6.2 The predicted variance 

n-Hh = /3(s-Ai s ) + e (62) 

(n- M/l ) 2 = /? 2 ( S -/i s ) 2 + 2/3( S - Ms )e + e 2 (63) 

E((n - fi h ) 2 \s) = (3 2 (s~ fi s ) 2 + v(s) (64) 

where v(s) = a + j3s and the expectation in the final step is again over realisations of the noise but 
not over realisations of the SST. The final equation is now slightly different than before because of the 
dependence of the variance on s. This has implications for the next step: 



J E((n- ^ h ) 2 \s)p{s)ds 


(65) 


J ((3 2 (s~^) 2 + v(s))p( S )d S 


(66) 


P 2 J (s~^ s ) 2 ds + J v(s)p(s)ds 


(67) 


(3 2 J(s- Li s ) 2 ds + J(a + (3s)p(s)ds 


(68) 


/3 2 <7 2 + (a + f3n s ) 


(69) 


P 2 v 2 s + Hh 


(70) 



We see that the part of the uncertainty in the hurricane number prediction that depends on the uncertainty 
in the SST-hurricane relationship now becomes /i^, the predicted mean number of hurricanes. In other 
words, the higher the mean number of hurricanes predicted, the greater the uncertainty (in absolute 
terms). 

Comparing the expression for the mean (equation l6T|) and the expression for the variance (equation l70|) 
we see that, in general, they are not the same. We conclude that the predicted hurricane distribution, 
although a mixture of poisson distributions, is not itself a poisson distribution. 

6.3 Standard errors 

The derivation for the standard errors is the same as for the linear-normal case: 

[i h = a + /3fi s (71) 

and so 

d/ih = da + fi s d(3 + /3dfj, s (72) 

and 

vax(fXh) = var(a) + /u 2 var(/3) + /3 2 var(/x s ) + 2/z s cov(a, 0) (73) 
As before, this is in fact exact. 

6.4 Linear sensitivity 

This is also the same as for the linear-normal case: 



(j, h = a + (3 [i s (74) 
dfi, 



P- - (75) 



and 



6.5 Summary 



Hh = a + /3n a (77) 
a 2 = (3 2 a 2 s +a + p f i s (78) 
var(^) = var(a) + /i^var(/3) + /? 2 var(/i s ) + 2/j, s cov(a, (3) (79) 

W. = " m 
w. = < 81 > 

Once again the fact that we assumed SST was normal was irrelevant... and again all we used about the 
SST was the information about the first two moments. 

6.6 Alternative representation 

The alternative representation is almost the same as for the linear-normal case: 

n = a + [3(s-s) + e (82) 

This then gives: 

fx h = a + 0(jji B - s) (83) 
The derivation for the standard errors is based on: 

dfih — da + (fj, s - s)d/3 + (3d^i s (84) 

giving: 

vax(^ h ) = var(a) + (p s - s) 2 var(/3) + /3 2 var(^ s ) + 2(^ s - s)cov(a, (3) (85) 
In this case the covariance term is not necessarily zero. 



7 The exponent ial-poisson model for the relationship between 
SST and hurricane numbers 

We now consider the third of our models for the relationship between SST and hurricane numbers, in 
which the mean number of hurricanes (given the SST) is given by an exponential of a linear function of 
SST, and the distribution of the number of hurricanes is poisson. 



n = exp(a + f3s) + e (86) 
n - Po(ratc = cxp(a + [3s)) (87) 

The variance depends on the SST, as for the linear-poisson model. 

There are now big differences in the following analysis because of the non-linearity. 

7.1 The predicted mean 

Because of the non-linearity, evaluating the mean number of hurricanes is now a bit harder. We actually 
have to do the integral, and so we actually do have to use the assumption that the SST distribution is 
normal, rather than just needing the first two moments as in the previous two models. 



/ 



exp(a + [3s)p(s)ds 



exp(a + [3s) 



2na. 



■cxp 



2a 2 



ds 



= exp( 



^ \f2Tta s I 



' ■■]> { [3s- {S 2 ^ s)2 ) 



(88) 
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= exp(a + f3(^ s + 0a 2 /2)) 
where we used the temporary variables: 
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2-ko s 
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ds 

ds 
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Ms + (3d 2 
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fi s -a 

0a 2 (2^ s ~0a 2 ) 



7.2 The predicted variance 

We can derive an expression for the variance in this case as follows: 

a 2 = E[(n-^ h ) 2 } 
= E(n 2 )-^ 2 

[i 2 we know already, so we just need to calculate E(n 2 ). 



En 2 = 



J E(n 2 \s)p(s)ds 
= J E[cxp(a + f3s) + e} 2 p(s)ds 
= J [cxp(a + f3s) 2 + v(s)]p(s)ds 
= j '[cxp(2a + 20s) + v(s)]p(s)ds 



= h+h 



h = J exp(2a + 2f3s)p(s)ds 

= exp(2a + 20(n + 20a 2 /2)) 

= [cxp( a + (3{ti + /?<7 2 /2))] 2 exp(/3 2 a 2 ) 

= /4exp(/3 2 a 2 ) 

h = J v(s)p(s)ds 

exp(a + 0s)p(s)ds 



= fJ-h 



so 



a 2 = E(n 2 )-fi (116) 

= fi 2 h cxp(/3 2 a 2 ) + fi h - n 2 h (117) 

= ^(exp(/3 2 <7 s 2 )-l)+ M (118) 

« HIP^+HH (119) 

where the final approximation will be accurate if /3 2 cr 2 << 1. 

7.3 Standard errors 

The derivation of expressions for standard errors follows the same logic as before, but is slightly more 
complicated because of the non-linearity: 
The mean is given by: 

M/l = e a+/J(^.+^/2) (120) 

= e x (121) 

where x = a + f3(p s + /3er 2 /2). 

dn h d^i h dnh, d^ h 
dHh = ~Q^ da + -Qp dl3 + ~d^~ d ^ s + ~&T » ( > 

= e x (da + (^+p t j 2 s )df3 + f3d^ s +l3 2 <j s d<j s ) (123) 

which gives: 

var(w.) = e 2a; [var(a) + ( Als +/3 ( T 2 )var(/3) + 2( Ms +/3a 2 )cov(a,/3) (124) 

+/3var(^i s ) + p 2 a s vav{a s ) + 2/3 3 a s cov(p s ,a s )} (125) 

Unlike the linear cases this is now not exact, and there is in fact an infinite series of higher order 
terms. Evaluating the second order terms would be a useful way to assess the accuracy of the linear 
approximation. 

7.4 Linear sensitivity 

We can also calculate the linear sensitivity as follows: 

= ^ (126) 
= /3e x (127) 

This makes sense: the sensitivity to the mean SST is proportional to the mean, and is proportional to j3. 
Also: 

^ = ^ (128) 
= 2 a s e x (129) 



or 

1 9 ^ h : p 2 e x (130) 



and we see that it is fractional changes in a s that matter. This is the only one of our three models where 
the standard deviation of the SST prediction a s has an impact on the mean hurricane numbers predicted. 
This arises because of the non-linearity: more uncertainty on the SST prediction leads to higher expected 
numbers of hurricanes. 



7.5 Summary 



H h = exp (a + /?(/*. +j9o*/2)) (131) 

a 2 h = M 2 (exp(/? 2 a 2 )-l)+^ (132) 

var(^) = e 22; [var(a) + (^ s +/3o- 2 )var(/3) + 2(^ s +/3fJ 2 )cov(a,/3) (133) 

+/3var(^i s ) + /3 2 a s vax{a s ) + 2/3 3 cr s cov(^ s , cr s )] (134) 



9cr. 



0e x (135) 
,8V,e* (136) 



7.6 Alternative representation 

The alternative representation is also slightly different than before: 

n = exp(a + /3(s - s)) + e (137) 

This then gives: 

/M =exp(a + /3(^-s + /?er 2 /2)) (138) 
The derivation for the standard errors is based on: 

dfi h = e x (da + (ii s -s + [3<T 2 )df3 + f3dn s +l3 2 (T s dcr s ) (139) 

giving: 

var(^h) = e 2x [var(a) + (jjl b - s + /?cr 2 )var(/3) + 2{fi s -s + /3a 2 s )cov(a, [3) (140) 

+/3var(^ s ) + (3 2 a s vax(a s ) + 2/3 3 a s cov([i s ,a s )} (141) 

In this case the cov(a, /?) term doesn't disappear, but will be much smaller than in the original formulation. 

8 The linear-poisson model for the relationship between basin 
and landfalling hurricane numbers 

We now consider the relationship between the number of hurricanes in the basin and the number at 
landfall. We will use our model for this relationship later when we consider the possibility of predicting 
landfalling hurricane numbers in a three step approach: by first predicting SST, then predicting basin 
hurricane numbers from SST, and finally predicting landfalling hurricane numbers from basin hurricane 
numbers. 

The only model we consider in this case is a linear-poisson model for the number of landfalling hurricanes 
as a function of the number of hurricanes in the basin: 

m = a' + (3'n b + e' (142) 
ni - Po(rate = a' + (3'n b ) (143) 

Comparing with the previous linear-poisson model described in section[6j the SST s has now become rib- 
There are only slight differences in the analysis compared with that model. 



8.1 The predicted mean 

Hi = E(n) =a' + (3' ii b 



(144) 



8.2 The predicted variance 



n - m = f3'(n h - (j, b ) + e' (145) 
(n-m) 2 = (3 ,2 {n b - n b f + 2(3'{n b - Mb )e + e' 2 (146) 
E((n- R )\n b ) 2 = (i' 2 (n b -ii b ) 2 +v(n b ) (147) 

(148) 



where v(n b ) = a' + fi'n b . 



a 2 = J2 E (( n -^ 2 \n b )p(n h ) (149) 

oo 

= ^(^K-^+^ntJJpK) (150) 



n=0 

oo 



= /3 /2 ^(n fe - M ) 2 + 5>(n 6 )p(n 6 ) (151) 

n=0 n=0 

= Z?' 2 ^ 2 + a' + $n b (152) 
= /3' 2 <7 6 2 +W (153) 

Once again, comparing the expression for the mean and the expression for the variance we see that, 
in general, they are not the same. The predicted hurricane distribution, being a mixture of poisson 
distributions, is not itself a poisson distribution. 

8.3 Standard errors 

fJH=a' + pi* (154) 

so 

dm = da' + fj, b d/3' + (3'd(j,b (155) 

and 

var(w) = var(a') + A^var(/3') + /?' 2 var(^) + 2fi b cov(a', f3') (156) 
As in the previous linear cases, this is exact. 

8.4 Linear sensitivity 

Hi=a' + Put (157) 

so 

dm 



dn b 



It doesn't make sense to consider |^ since a 2 = fi b . 
8.5 Summary 



Ht = a' + l3'fi b (159) 

a 2 = P' 2 <T 2 b +m (160) 

var(^) = var(a')+M 2 var(/3') + /3' 2 var( Mb ) + 2 J u fc cov(a',/3') (161) 

|^ = (162) 



8.6 Alternative representation 



ni = a' + f3'{n b -n£) + e' (163) 
ni ~ Po(rate = a' + /3'(rib - (164) 

(Ml = a' + f3'(nb -nb) (165) 
(f/ij = da' + (/if, — ri b )dj3' + P'dfj, b (166) 
var(/x ; ) = var(a') + (fi b - n£) 2 var(,S') + /3' 2 var(^ fc ) + 2{fi b - n£)cov(a', /?') (167) 

9 Predicting landfalls from basin numbers, and basin numbers 
from SST 

The models given in sections [5j [6] and [7] for relating SST to hurricane numbers (which we now take as the 
number of hurricanes in the basin) can now be combined with the model given in section |8] that relates 
the numbers of hurricanes in the basin to the number at landfall. 

We consider two cases below, based on the linear-poisson and exponential-poisson models for the number 
of hurricanes in the basin, and combined with the linear-poisson model for the number of hurricanes at 
landfall. 

9.1 Linear poisson model from SST to basin, linear-poisson model from basin 
to landfall 

9.1.1 Mean and variance 



From section [B~5l we see that given an SST distribution N(/i s , a 2 ) the mean and variance of the distribution 
of the number of basin hurricanes is: 

[i b = a + I3(i s (168) 

o\ = (3 2 <j 2 + t i b (169) 

and from section [8751 we see that given the mean and variance of the number of basin hurricanes the mean 
and variance of the number of landf ailing hurricanes is: 

Hi = a' + 0'fi b (170) 

af = /3' 2 a 2 + M; (171) 

Putting these together in order to get from SST to landfalls in one step, we get: 

Hi = a' + p'nb (172) 

= ct+p(a + Pn.) (173) 

= a 1 + /3'a + P'0n s (174) 

a? = P' 2 a 2 b +m (175) 

= /3' 2 (/3 2 a 2 +^)+/i, (176) 

= /3' 2 /3 2 a 2 + (3' 2 n b + m (177) 

9.1.2 Standard errors 

For the SST to basin hurricanes part we have: 

var(^i fc ) = var(a) + ^ 2 var(/3) + /3 2 var(/i s ) + 2^ s cov(a, (3) (178) 
and for the basin to landfall part we have: 

var(/ij) = var(a') + ^var(/3') + /3' 2 var(/i 6 ) + 2^ fc cov(a', ff) (179) 
and these two expressions can easily be combined to give a one-step expression for var(///). 



9.2 Exponential poisson model from SST to basin, linear-poisson model from 
basin to landfall 

9.2.1 Mean and variance 

From section [7. 5 1 we see that given SST distribution N(p s , a 2 ) the mean and variance of the distribution 
of the number of basin hurricanes is: 

fi b = exp(a + /3(^+/3a s 2 / 2 )) (180) 
o\ = ^(exp(/3 2 ( r s 2 )-l) +Aib (181) 

and from section 18.51 we again see that given the mean and variance of the number of basin hurricanes 
the mean and variance of the number of landfalling hurricanes is: 

111 = a' + p'm (182) 
a? = I3 l2 a 2 b +m (183) 

Putting these together: 

IH = a 1 + /?'exp(a + /3( Ms + /?cx 2 /2)) (184) 
a? = /3'V 2 (exp(/3 2 a 2 )-l)+/3'V, l + W (185) 

9.2.2 Standard errors 

For the SST to basin hurricanes part we have: 

var(/i&) = e 2x [va,r(a) + (ji + /3<r 2 )var(/3) + /3var(^i) + /3 2 crvar(cr) + 2{(i + /3a 2 )cov(a, 0)] (186) 
and for the basin to landfall part we have: 

var(/x;) = var(a' ) + n 2 b v&r((3') + /3' 2 var(^ b ) + 2(i b cov(a', (3 f ) (187) 
and again these two expressions can easily be combined to give an expression for var(/^). 
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